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Semiconductor crystals such as Hgi-«CdxTe grown by unidirectional solidification 
Bridgmann method have shown compositional segregations in both the axial and radial 
directions (Lehoczky et. al. 1980, 1981, 1983). Due to the wide separation between the 
liquidus and the solidus of its pseudobinary phase diagram (Lehoczky and Szofran 1981 ), 
there is a diffusion layer of higher HgTe content built up in die melt near the melt-solid 
interface which gives a solute concentration gradient in the axial direction. Because of the 
higher thermal conductivity in the melt than that in the crystal there is a thermal leakage through 
the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the 
radial direction. Hart (1971), Thorpe, Hutt and Soulsby (1969) have shown that under such 
condition a fluid will become convectively unstable as a result of different diffusivities ot 
temperature and solute. It is quite important to understand the effects of this thermosolute 
convection on the compositional segregation in the unidirectionally solidified crystals. To reach 
this goal, we start with a simplified problem. We study the nature of fluid flows of a stratified 
solution in a cylindrical container with a radial temperature gradient. The cylindrical container 
wall is considered to be maintained at a higher temperature than that at the center ot the 
solution and the solution in the lower gravitational direction has higher solute concentration 
which decrease linearly to a lower concentration and then remain constant to the top of the 
solution (see Fig. 1). The sample solution is taken to be salt water with its properties listed m 

table 1. 


Table 1. Material Properties for 3.5% salt water 


Mass diffusivity for salt water 
Thermal diffusivity 
Density of salt water 
Viscosity of salt water 
Kinematic viscosity 
Prandtl number 
Schmidt number 
Volume xpansion coefficient 

( thermal) 

Diameter of the cylinder 
Grashof number for temperature 


«c = 1 .3E-9 m 2 /s 
cxt = 1 .5E-7 m 2 /s 

P = 920 K g/ m3 
|i = 7.6E-4 Kg/m-s 
t> =8.26E-7 m 2 /s 
P^tVoer =5.65 
Sc=fVtXc = 635.4 

fc =2.75E-4 
R = 0.005 m 
Grp = 22.23 A 2 


Under the Boussinesq approximation, the vector form of the governing equations for 
incompressible fluid flow of the system are: 


Momentum Conservation 


^ + mVu = -IVp - g{ WT-To) +Pc (C-Q> ) Fz + uV 2 u ^ 

Continuity Equation 

V*u = 0 [2] 

Energy Conservation 

f + 5vr-^ (3] 
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Solute Conservation 


^ + iT»VC = otcV 2 C 


[ 4 ] 


Where u is fluid flow velocity, T and C are the temperature and the solute concentration of the 
fluid, respectively, and a T and are the diffusivity for temperature and solute, respectively. 

For axial symmetric boundary conditions, the governing equations in cylindrical 
coordinates are: 
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§ +u § +l0 |c (*c +i ^c + *c 

3t 3r 3z I a r 2 r 3r dz 2 
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The equations can be nondimensionalized by setting the temperature and solute 
concentration as 


T =T 0 + ATT* 


and 


C =C 0 + ACC*, and 


scaling the va riables by a factor F; i. e V = F V* Scaling length by R, velocity by 

YgprATR, time by R/fgprM, pressure by pgPrATC temperature by AT and solute by AC. 
After the scaling and dropping all the *, the dimensionless equations become: 
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[ 10 ] 


Grr'^+u^ + +1 ^ 

\ 3t 3r dz I dz 1 3r 2 r 3r dz 2 


+ Gr T 1/2 ( T + (Grc/Grj) C) , 


[ 11 ] 
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IdT 5T 3T\ d 2 T i 3T 

PrGrTl/ 1aT +u dF +CD aTr ar 2 r * dz 2 , 
m/ac ac aq 32 C + i^£ + §^. 

ScGrr 1/2 |^ + u ^ + 00-^)= -^ + f * + ^ 


Where the thermal and solutal Grashof numbers respectively , Git, Grc are defined by: 


GlT jfe*™i 


Gr c =Sfc^i 


v 2 and u 2 


The Prandtl number, Pr, and Schmidt number. Sc, are defined by: 


Sc=tVac and Pr=War. 

These governing equations show that the flow characteristic are determined 
uniquely by Git, Grc, Pr and Sc. These equations have been solved by the FIDAP program 
developed by Fluid Dynamics International, Inc. Preliminary results reveal that any convection 
influences the solute field much more than it does the temperature field. Since the Prandtl and 
Schmidt numbers are properties of the particular fluid, the variations with Git and Grc were 
primarily studied. Figure 1 through 4 depict the computational domain and some of the results 
obtained. 



Fig. 2. Concentration and temperature fields for 
Git =0.04 and Grc/Gir =1.0 
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Fig. 3. Concentration and temperature fields for 
fields for . Grp =0.49 and Grc/Grp =1.0 


Fig. 4. Concentration and temperature 
Grr =0.49 and Grc/Grp =1000 
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